Ministere de I'education de la Federation de Russie 
Universite d'etat de I'Oural 

Faculte de mathematiques et mecanique 
Chair d'analyse mathematique et theorie des fonctions 



LA MEILLEURE APPROXIMATION MOINDRES CARRES 
PAR LES L-SPLINES DE DEUXIEME ORDRE 



These pour I'obtention du magistere 
speciahte « Mathematiques, mathematiques apphquees » 



Chef de la chair presentee et soutenue par 

professeur Arestov V.V. Kaioumov Alexandre Rachidovitch 



Directeur de these — 

membre correspondant de I'Academie des Science de la Russie, 

professeur, docteur es sciences physico-mathematiques, 

chef du Departement d' approximation 

de rinstitut de mathematiques et mecanique 

de I'Academie des science de la Russie 

Soubbotine Youri Nikolaievitch 



Ekaterinebourg, 2003 



Resume 



Mots-clefs : L-SPLINE, SPLINE EXPONENTIELLE, SPLINE TRIGONOMETRIQUE, MEILLEURE 
APPROXIMATION MOINDRES CARRES, PROJECTEUR ORTHOGONAL, CONSTANTE DE LE- 
BESGUE. 

Sujet de la these. Dans cette these nous etudions I'operateur de Tapproximation moindres carres 
par Ics i-splinos gonoroos par im oporatour difFerentiel linoairo do douxiomc ordro — autromont dit, 
I'operateur de projection orthogonalc sur I'cspace de ces L-spUnes. Nous considerons coninie I'operateur 
differential les operateurs generant les L-splines exponentielles et trigonometriques. Nous traitons la 
question suivante : la norme uniforme d'un tel operateur de projection, est-elle bornee uniformement en 
tous les partitionnements du segment ? 

Resultats principaux. Nous demontrons que les operateurs de projection orthogonale sur de diffe- 
rents espaces de i-splines sont bornes uniformement en tous les partitionnements du segment. Pour les 
spHnes exponentielles avec les exposants d'une m6me valeur absolue, une estimation exacte est obtenue 
pour la norme du projecteur. Pour les splines exponentielles avec les exposants arbitraires de signes diffe- 
rcnts, nous demontrons que le projecteur est uniformement borne. Pour les splines trigonometriques avec 
des contraintes habituelles sur le segment, une estimation est obtenue pour la norme du projecteur. Nous 
traitons aussi un exemple modele dans le but de montrer que ces contraintes sur le segment ne sont pas 
du tout necessaires pour que le projecteur orthogonal soit borne. 

Methodes employees. Nous examinons les matrices de Gram des bases de i-B-splines, aussi que 
leurs inverses. Dans le cas des spHnes trigonometriques, au cours de I'etude de I'exemple modele, nous 
etablissons un phenomene d'interference, qui permet a I'approximant trigonometrique de rester borne, 
tandis que ses coefficients et les splines de base tondont vers I'infini. 

Histoire du probleme et portee des resultats. L'etude des normes des projecteurs orthogonaux 
sur de differents espaces fonctionnels est un probleme classique de la theorie d'approximation bien im- 
portant on lui-m6me. Applique aux splines polynomials, ce probleme est connu sous le nom de probleme 
de Carl de Boor. 

Dans le cas des splines, la norme du projecteur orthogonal sur I'espace des splines de degre k — 1 est 
liee a la norme de I'operateur d'interpolation par les spHnes de degre 2k — 1, ce qui accentue I'importance 
de son etude. Professeur C. de Boor a lance une etude systematique de ce probleme, qui I'a amene a 
la conclusion que la norme du projecteur devait etre bornee independamment du partitionnement du 
segment, une conjecture qu'il a formulee en 1972. 

Malgre de nombreux resultats partiels, le probleme de C. de Boor restait ouvert pendant plus de 25 
ans, jusqu'a ce qu'il n'ait ete resolu positivement par Dr. A.Yu. Chadrine en 1999. 

Le probleme considere dans cette these constitut done la premiere etape dans la generahsation du 
probleme de C. de Boor aux L-spHnes. 



1 Introduction 



1.1 Position du probleme 

Nous introduisons les notations suivantes. Fixons un segment [a, b] et designons par 

- A = {a = ti < t2 < ■ ■ ■ < t„-i < tn = b} — un partitionnement du segment [a, 5], 

- Ai = [ti,ti+i] (i = 1, . . . , n — 1) — les sous-segments de ce partitionnement, 

- L = D'' + Gk-iD''^^ + • ■ • + aiD + qqI, a; e R — un operateur differential lineaire d'ordre k a 
coefficients constants, 

- Nl = {/ e C'^[a, 6] : Lf{x) = 0, a; G [a, 6] } — le noyau de I'operateur L. 

Alors I'espace des L-spUnes genere par le partitionnement A et I'operateur differentiel L est defini par 

S(i,A) = {s e C''-^[a,b] : s|a. e Nl, z = 1, . . . , ?i - 1 } . 

Considerons I'operateur P§ = Ps(l.a) de projection orthogonale (relative au produit interieur dans 
I'espace L2[a,b]) de I'espace C[a, 6] sur son sous-espace S (L, A) : 

Ps:C[a,5]9/^Ps/e§(L,A), 

ou Pg / est uniquement defini par I'une quelconque des deux conditions suivantes : 

1. (/,s) = (Ps/,s), V,seS(L,A), 

2. II/-PS/II2- min 

sGa (L, A) 

Nous nous occupons de la norme de I'operateur Pg a titre d'un operateur de C[a, b] a valeurs dans C[a, b] : 

Hi's Hoc = sup { \\Psf\L:feC[a,bl \\f\\^^l}. 

II est aise de voir que lorsque le partitionnement A est fixe, I'operateur Ps(l,a) est borne, c'est-a-dire a 
une norme finie. Nous sommes done naturellement amene a la question suivante : la famille d'operateurs 
de projection Ps(l.a), est-elle bornee uniformement en tons les partitionnements A du segment [a,b] ? 



1.2 Histoire du probleme 

Le probleme que nous venons de formuler est la generalisation au cas des L-splines du probleme 
correspondant pour les splines polynomials : on obtient le dernier en posant L = D'^ dans le premier 
(les L-splines ainsi obtenues comcideraient avec les splines polynomials d'ordre k, de degre fc — 1). 

Faisons un expose de I'histoire du probleme pour les splines polynomiales. L'hypothese que les normes 
uniformes des projecteurs orthogonaux Ps^[\) sur les espaces §/c(A) des splines polynomiales sont bornees 
independamment des partitionnements A du segment [a, b] a ete enoncee par Carl de Boor dans en 
1972. C'est pourquoi cette liypothese etait connue sous le nom de probleme de C. de Boor. Designons 

Cfe =SUp||Ps,(A) 
A 

Dans le cas de fc = 2 (splines lineaires ou lignes brisees), Z. Ciezielsky 0j a montre en 1963 I'estimation 
C2 ^ 3. K.I. Oskolkov a obtenu en 1979 un resultat d'ou il decoule que C2 = 3. 

Dans le cas de fc = 3 (splines paraboliques) C. de Boor ^ a montre en 1968 I'estimation C3 < 30. 

Dans le cas de fc = 4 (splines cubiques) C. de Boor a annonce en 1979, en partant de resultats 
d'experiences numeriques, I'assertion que C4 ^ 2|5 ~ Autant que je sache, une demonstration exacte 
de cette assertion n'a jamais ete publiee. 

Apres cela, il n'y avait eu aucun progres pendant 20 ans — jusqu'au 1999, quand le probleme de 
C. de Boor a ete resolu par Dr. A.Yu. Chadrine dans i2i- H a verifie que Ck < +00 (fc > 2), bien qu'il 
n'est pas arrive a obtenir une estimation superieure constructive. En utilisant la methode d'Oskolkov, 
il a neanmoins obtenu une estimation inferieure : Cfc ^ 2fc — 1. En partant de resultats d'experiences 
numeriques et de quelques considerations theoriques, il a enonce la conjecture que Cfc = 2fc — 1 (fc ^ 2). 

Autant que je sache, le probleme correspondant pour les L-splines n'a jamais ete examine. 
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1.3 Formulation des resultats 



Dans cette these nous nous limitons a I'analyse de L-splines de deuxieme ordre de type special, a 
savoir des splines exponentielles et trigonometriques de deuxieme ordre. Quoique celles splines-ci soient 
analogiques aux splines polynomiales lineaires (lignes brisees) qui ne presentent aucune difEculte, le 
passage des splines polynomiales aux L-splines rend la situation assez complexe. 

Les classes suivantes de L-splines sont traitees et les theoremes suivants sont demontres. 

1. Splines exponentielles avec les exposants d'une meme valeur ahsolue. Considerons les L-splines 
generees par I'operateur differentiel L = — A^. Une telle spline coincide sur chaque sous-segment 
du partitionnement avec une combinaison lineaire des fonctions e^^'-^ et e"'*'^. 

Theoreme 1. Pour les L-splines generees par I'operateur differentiel L = — on a : 

c(+A, -A) sup II Ps(L,A) lloo ^ 3- 

Notons que la methode employee pour montrer I'exactitude de I'estimation ci-dessus nous fournit 
une autre demonstration du resultat d'Oskolkov (que I'estimation par Ciezielsky pour les splines 
polynomiales lineaires est exacte). 

2. Splines exponentielles avec les exposants de signes dijferents et de differentes valeurs absolues. Consi- 
derons les L-splines generees par I'operateur differentiel L — {D — Xi){D — A2), ou A2 < 0, Ai > 0. 
Une telle spline coincide sur chaque sous-segment du partitionnement avec une combinaison lineaire 
des fonctions e^''' et e^"^. 

Theoreme 2. Pour les L-splines generees par I'operateur differentiel L = [D — \2){D — ou 
A2 < 0, Ai > 0, on a : il existe une constante C(Ai,A2) qui ne depend qu'au nombres Ai,A2 et ne 
depend pas an segment [a, b], telle que la famille d'operateur de projection Ps{l.a) est bornee par la 
constante C(Ai, A2) uniformement en tous les partitionnements A du segment [a,b] : 

c(Ai,A2,a,5) =^sup||Ps(£^A) |L C'(Ai,A2). 



3. Splines trigonometriques. Considerons les L-splines generees par I'operateur differentiel L = D^+X^. 
Une telle spline coincide sur chaque sous-segment du partitionnement avec une combinaison lineaire 
des fonctions sin Ax et cos Ax. 

Ce cas-ci differe des deux cas precedents et du cas des spHnes polynomiales : comme le probleme 
d'interpolation lagrangienne par les fonctions trigonometriques devient insoluble sur un segment 
d'une longueur de tt/A, les proprietes des L-splines de base empirent lorsque le diametre du par- 
titionnement tend a cette valeur. C'est pourquoi en traitant les splines trigonometriques on essaie 
habituellement d'eviter la proximite du diametre du partitionnement a la valeur tt/A. De plus, on 
introduit souvent une contrainte encore plus forte : d{A) < ^ (ou d{A) — max ^{ti+i — ti) — 

est le diametre du partitionnement A). (Voir, par exemple, jH], page 455, oil 0, page 135.) 
Designons 

Af(A,T)=^ sup ||Ps(L,A)|L- 
A-d(A)<T 



Theoreme 3. Pour les L-splines generees par I'operateur differentiel L = + A^ on a : 
M (A, r) «C 



(1 — cosr)^ „ TT 

; siO<T^-, 

(T — smrj smr 2 

(1 — cos r)^ TT 

si — < T < n. 



(t — sin r) sin r 2 
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Mais quand r 



vr — le deuxieme membre de I'estimation ci-dessus est equivalent a 



8 1 



TT (tT — t)2 



Une question se pose d'une fagon naturelle : est-ce que cette asymptotique correspond a la conduite 
actuelle de la valeur M (A, r) ? En general cette question semble difficile, et nous considerons ici un 
exemple modele, quand parmi tous les sous-segments il n'y a qu'un d'une longueur proclie a j. 

Theoreme 4. Soit [a,b] = [0, un des sous-segments ait une longueur proche d j, la somme des 
longueurs de tous les autres sous-segments soit e. Alors 



Les resultats ci-dessus ont ete publics dans |!?l llfll[Tnil2j . 
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Theoreme 5. Soit [a,b] = [0, f]. Alors 



M (A, tt) < +00 
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